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1. Introduction 

At the turn of the century the study of the quantum manifestations of classical chaotic 
systems suffered a significant change. Before, the spectral statistics were amply discussed 
and it was shown that they follow the Random Matrix Theory (RMT) predictions |1|, 0]. 
The study of wave function properties, however, presents inherent difficulties arising 
from the dependence on the basis used, forcing either to specify one or to define basis 
independent quantities. Recent progress has been made in this respect in the study of 
average properties of eigenstates 0, || and of the significant statistical deviations 
from RMT [P, |^. Nevertheless, not much systematic work exists on the eigenfunction 
fiuctuations in dynamical systems 0. In this work we shall contribute on the last 
subject. 

In a recent paper |^ it was found that the suitably averaged matrix elements 
between the eigenfunctions (EF) of two arbitrary Hamiltonians Hq and H are well 
described in the semi-classical regime by a classical phase-space integral. Specifically, 
if we define (pa, to be the eigenfunctions and eigenvalues, respectively, of i^o and 
ipi,Ei those of H, one finds to a good approximation 



where g{e,E) is given by 

gie,E) = ldpdqS{Ho{p,q)-e)SiH{p,q)-E), (2) 

which was called in reference the classical eigenfunction for fixed E. Here po{e) is the 
level density of Hq calculated by means of Weyl's formula. By the symmetry of equation 
(|^), the local density of states(LDOS) can be calculated using the energy density of H 
instead of Po{e) and maintaining fixed e in equation (0). For the details, in particular 
about the way in which the l.h.s. of equation (|1|) must be averaged to obtain meaningful 
results, see reference f^. This study was exemplified by two systems of anharmonic 
oscillators in one dimension, one coupled and the other uncoupled. In the present letter, 
our interest is focused on the fiuctuations of the quantum-mechanical wave functions 
around this classical limit in the chaotic case. 

2. A Random-matrix model 

In the usual description of chaotic systems by random matrices, the restrictions implied 
by equations (|l]) and (@) are not present. Rather, one attempts to deduce the average 
properties of the eigenfunctions given the structure of the random matrix ensemble in 
some particular basis. Instead, here we choose pairs of matrices {H, Hq) of size N x N 
in such a way that the condition 

{m^p,)f) = u, (3) 

is always fulfilled, but the matrices are otherwise arbitrary. The angular brackets denote 
the average over the ensemble of matrix pairs and lij stand for any numbers given by 
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outside constraints such as Under these circumstances, we wish to determine the full 
distribution of the matrix elements {(pilipj). We then proceed to compare the predictions 
of this random matrix model with numerical results on models similar to that studied 
in 0. 

To solve the above problem, we let ourselves be guided by the following 
considerations: the quantities we need to model, namely the {(pilipj), are nothing else 
than the matrix elements of an orthogonal matrix (a unitary one in the case where time- 
reversal invariance is broken, but this does not affect our conclusions). We therefore 
need a random matrix model for orthogonal matrices with prescribed expectation values 
for the intensities lij. Note that, in the large limit, the Haar measure over the group 
of orthogonal matrices can be replaced, up to corrections of order by independent 
Gaussian distributions for all matrix elements Oij, all having a variance In 
other words, what we need is a random matrix model where the average intensity 
lij = is given. If Jj j <^ 1, we can consider as an effective dimension and 

we expect to get the desired result, up to corrections of order 0{lij), if we replace the 
average by a simple Gaussian average. We may then, to this level of accuracy, take the 
Oij as independent Gaussian variables with variances given by lij. 

Note that we postulate a distribution for orthogonal matrices with the correct values 
for lij. We do not actually derive this distribution, but simply verify that it has all 
the required properties. If the classical wave function takes very large values and the 
eigenfunctions of H do not have a sufficiently large number of components in the basis 
of Hq, it may happen that some lij ~ 1. Clearly, for these matrix elements the model 
will not apply; we shall see later that this happens near peaks or singularities in the 
classical wave function, but then we cannot really make any comparison with the specific 
system anyway. However, for the vast majority of matrix elements, we can expect that 
the amplitudes are Gaussian distributed and if we divide the amplitude Oij by ^Jhj we 
will find a standard Gaussian. 

What deviations from the above predictions should we expect from a theoretical 
point of view? Clearly, scars produce an excess of very small amplitudes, because 
a few exceptionally large amplitudes pick up more of the intensity than expected from 
the classical calculation. Would we also see the large amplitudes? Probably not in 
a statistical analysis against our model, because these will mainly occur in the region 
where the classical function is large and we will usually exclude this region: the condition 
lij <^ 1 is violated there, unless we reach very high spectral densities, which is scarcely 
possible in a numerical experiment. In a real experiment, resolution might well make 
such a high-density region inaccessible also. If, on the other hand, localization occurs 
due to disorder or due to the fact that the system does not cover the whole phase space 
on the Heisenberg time scale then we may indeed also see irregularities beyond the 
realm of very small amplitudes. 

However, all deviations mentioned above should only be important if one of the 
Hamiltonians, say Hq, is integrable in the classical limit. If both are chaotic, and we 
exclude situations in which the two Hamiltonians are, in some sense, closely related, we 
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can expect not to see any effect of the scars in the amphtudes. The reason for this can 
be understood in terms of the traditional picture due to Berry of the eigenf unctions 
in phase space: For a chaotic Hamiltonian, they are expected to cover phase space 
essentially in a uniform way, up to rather small concentrations on periodic orbits. In 
integrable systems, on the other hand, eigenfunctions are localized on well-defined tori, 
with only half the dimension of the full phase space. The overlap between two chaotic 
states is therefore far less likely to become anomalously large than the one between an 
integrable state and a chaotic one. 



3. Numerical results 



We now test the Gaussian property against anharmonic oscillator models. We shall 
choose the expansion of a chaotic system in terms of an integrable one; in particular, 
we have chosen two particles in a quartic oscillator potential. This ensures a system 
with scaling properties, for which the classical properties do not change as a function 
of energy. We restrict our attention to antisymmetric wave functions since in this case 
we reach the semi-classical limit much more rapidly than for the symmetric case. The 
calculation is performed using the basis of the uncoupled oscillators, which in turn we 
approximate in a harmonic oscillator basis The Hamiltonian used is 

n 2 n n n 

^ = Ey + «E^' + /^ E + 7 E h4 + ^'^i] ' (4) 

i=l ^ 1=1 l<i<j l<i<j 

where in this case n is equal to two. We have also considered the case n = 4 with overall 
similar results |]10|. We shall use two Hamiltonians: one with the same parameters as 
in [Q, namely a = 10, (3 = —5.5 and 7 = 5.6, which we call Hi, the other with the 
parameters a = 10, /? = 7 = —4.15, which we call H2. The Hq Hamiltonian will have 
a = 10 and /? = 7 = 0. 

We analyze the eigenfunctions in terms of the classical eigenfunction as given in 
equation ( ^ at fixed energy E. The integral is calculated by the Monte Carlo method. 
We find very good agreement as shown in figures and |[ There the classical EF and 
an average over 101 EF's of the perturbed Hamiltonian are plotted using the method 
of reference Note that the quantum functions are not reliable at the upper end of 
the classical energy range of Hq for the higher lying states, although their energies are 
quite reliable. At the lower end of the spectra, on the other hand, the amplitudes are 
very good, and we find a consistent approximation to an exponential decay of intensities 
in the classically forbidden region as shown in figure |l], with some system dependent 
oscillations (these disappear in the 4-particle case pO|] ). 

We now proceed to analyze the amplitude fluctuations. We do this in the wings 
of the wave functions far from the peak, in regions where the classical function varies 
slowly and is sufficiently small to ensure lij <^ 1; of course we restrict our attention 
to reliable amplitudes. To this end, we first cut out the parts of the wave function 
which are either too high in energy so that they are not reliable, or which lie outside 
the classically allowed region. We further cut 4 states on either side of the singularity 
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at the peak of the classical EF. We do this because the fluctuations around the peak 
are large and the peak itself at energy ii^ is a singularity of the classical EF. We set the 
norm of the rest to one in both EF and the classical EF. Then we proceed to unfold the 
EF dividing the quantum EF by the classical one, deflned in equation ( |1]). In order to 
compare the different unfolded EEs we renormalize them again. 

The corresponding regions in the wings, labeled B and C in flgure ^ are the ones 
with the best quantum-classical correspondence. We use the intensity shape instead of 
amplitudes for clarity, but all the calculations were performed on the latter. In order to 
avoid the rapidly varying region, for the cases shown below we drop a window of 4 mean 
energy level spacings A centered in the eigenfunction for Hi and one of 80 A for H2] the 
end of the C region is 400A away from the center for both Hamiltonians (The values 
of A are, respectively, of 0.826 and 0.729.) As we cannot perform ensemble averages, 
we will perform energy averages within these windows after dividing the amplitudes by 
the square root of the local average intensity obtained from the classical function, which 
agrees well with the quantum average. As the center of each EF changes in energy, 
the window center changes but its width remains constant. The amplitude distribution 
we flnd, is plotted in flgure § for the superposition of the results of regions B and C 
on 101 EEs; for low- lying states the shape is far from Gaussian while for high- lying 
states we flnd fair agreement with the Gaussian behavior except for the excess of small 
intensities, which we expect due to scars. Such scars were seen in the two-body system as 
exceptional states with much narrower intensity distributions and smaller participation 
ratios similar results are found for the four-body system |]T0[. Nevertheless, a 
semi-log plot of the amplitude distribution shows a good parabolic shape in the wings, 
even for the zone C in a low- lying states around ipiQo{see flgure 0(b)). 

We now test our assumption that scar effects are not seen when we expand the 
chaotic Hamiltonian in a basis of another chaotic Hamiltonian, instead of an integrable 
one. For such an expansion we put the EF of the H2 Hamiltonian in terms of the Hi 
EEs. The quantum-classical correspondence is shown in flgure In this case the 
exponential decay in the classically forbidden area shows a hump, for which we have 
no explanation. The B zone is wider and in consequence the statistics are better as we 
show below. The amplitude distribution in the same region as in the previous case flts 
the Gaussian better, as shown in flgure ^. The excess of small amplitudes decreases and 
the agreement is better in a wider energy regime. A similar result is observed if we drop 
localized states in the statistics for the previous case. Beyond all these features, if we 
consider small windows in the tail of eigenfunctions we flnd statistically good Gaussians 
for both cases. In figure |^ we show some of them. The window width is of 20 mean level 
spacings in order to have a sufficient number of amplitudes (~ 1000) of the 101 EFs 
considered for the average. They have energies between 1000 and 1020 in figure ^(a) 
for the state 900 of Hi and from 640 to 660 for the state 500 of H2 in figure |^(b). The 
fluctuations in figure ^ are larger than in the previous figures, but all of them are inside 
the statistical deviation, as shown by using the test per bin, which is Xw = 12.289 
for (a) and Xw ~ 14.483 for (b). For clarity we plot the histograms with larger bins and 
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normalized to (total number of events) x (bin width). We cannot get such a good fit to 
the Gaussian distribution for all energy ranges; the larger the window width, the worse 
the observed fit. 

4. Conclusion 

We have analyzed the fluctuations of quantum-mechanical eigenfunctions with respect 
to their classical limit. Using a simple random-matrix model, the amplitudes are shown 
to follow a Gaussian distribution. This is confirmed by a numerical calculation using 
systems of two particles interacting through anharmonic potentials; agreement improves 
as we move up in the spectrum. We further find evidence for scars in an excess of small 
amplitude values as compared to the theoretical prediction if we express the cigcnstate 
of the chaotic Hamiltonian in term of an integrable one. This effect decreases markedly 
when both Hamiltonians have chaotic dynamics. 
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Figure 1. Averaged eigenfunction around (a)V'5oo and (b)V'9oo both for the Hi 
Hamihonian with the parameters given in the text, with -E500 = 579.267 and 
-£'900 = 846.680, respectively. The smooth curve represents the corresponding classical 
EF. Note the general exponential decay plus oscillations in the classically forbidden 
region. 




Figure 2. Eigenfunction tail zones B and C considered for analysis. The extreme 
cutoffs are determined by the quantum classical agreement, and renormalized for 
unfolding. The center window is 4A and 80A width for the corresponding cases 
analyzed here. 
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Figure 3. Amplitude distribution for the unfolded EF with energies -Eioo = 204.013 
(a)-(b) and Sgoo = 846.680 (c)-(d) in zones B and C. The Gaussian with the same 
moments as the distributions are indicated by the continuous line. Note that (b) is 
plotted in semilogarithmic scale. 




Figure 4. Averaged eigenfunction around -0500 for H2 Hamiltonian in the Hi basis 
with the parameters given in the text and at energy E^qq = 611.717. The smooth 
curve represents the corresponding classical EF. 
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Figure 5. Amplitude distribution for the unfolded EF £^500 = 611.717 of H2 
Hamiltonian in the zone B (a) and C (d). The Gaussian with the same moments 
as the distributions are showed by the continuous line. 
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Figure 6. Histogram of the amplitudes distribution in a window of width of 20 
energy units for 100 EFs around (a)'i/'9oo of Hi and (b)'i/'5oo of H2. The normalization 
corresponds to the number of events per bin, the total number of events being 2525 
and 2595 for (a) and (b), respectively. The corresponding values of are Xie = 12.289 
and X16 = 14.483. The continuous line corresponds to a Gaussian curve with the same 
moments as the histograms. 



